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Abstract 

Using the chiral representation for spinors we present a particularly transparent way to generate the 
most general spinor dynamics in a theory where gravity is ruled by the Einstein-Cartan-Holst action. In 
such theories torsion need not vanish, but it can be re-interpreted as a 4-fermion self-interaction within 
a torsion-free theory. The self-interaction may or may not break parity invariance, and may contribute 
positively or negatively to the energy density, depending on the couplings considered. We then examine 
cosmological models ruled by a spinorial field within this theory. We find that while there are cases for 
which no significant cosmological novelties emerge, the self-interaction can also turn a mass potential 
into an upside-down Mexican hat potential. Then, as a general rule, the model leads to cosmologies 
with a bounce, for which there is a maximal energy density, and where the cosmic singularity has been 
removed. These solutions are stable, and range from the very simple to the very complex. 

1 Introduction 

The greatest tragedy of XX century physics was by and large that gravity refused to partake in the successes 
of quantum field theory and the gauge principle. This has led to numerous schemes purporting to supersede 
both classical relativity and standard quantization, but none of them was fully successful. A very basic 
question can be asked: if gravity is to be seen as a gauge theory, which symmetry group is being gauged? It 
is possible to regard General Relativity as the symmetry broken phase of a gauge theory of groups for which 
the Lorentz group is a sub-group. Particular examples are the Poincare group [j][2] and the de Sitter/anti-de 
Sitter groups HHG3. 

A commonality of all these approaches is that they recover the so-called 'spin-connection' formulation of 
General Relativity, in which the gravitational field is described by two independent ingredients. The first, 
referred to as the spin-connection, acts a gauge field for the Lorentz group whilst the second is a Lorentz 
valued spacetime one-form called the co-tetrad. It is from the latter that the familiar metric tensor may be 
constructed. This formulation is elegant and desirable, but it does open up the doors to space-time torsion. 
In the presence of spinors one is naturally led to build actions directly dependent on the spin-connection. 
These produce a source term in the torsion equation of motion: and thus torsion is forced upon gauge 
theories of gravity, whenever spinors are present. 

It turns out that, at least in the minimal theories with this feature, the torsion is algebraically related to 
the spin density. Therefore it can explicitly be integrated out of the theory, at least classically. It is found 
that the theory is equivalent to a torsion-free theory endowed with a 4-fermion self-interaction. However, 
even without considering more elaborate theories, with propagating torsion for example, one has to face a 
number of different possibilities in the fermionic couplings. Terms which usually are boundary terms no 
longer drop out of the equations whenever torsion is present. Therefore these theories present a richness of 
possibilities, and the question naturally arises as to how to constrain them. 

The first part of this paper (Sections [5] and deals with the formal aspects of this matter. Availing 
ourselves of the chiral representation for spinors, we present a particularly transparent way to generate the 
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most general spinorial dynamics in a theory where gravity is ruled by the Einstein-Cartan-Holst action. 
We also work out explicitly the 4-fermion self-interaction in the equivalent torsion-free theory. The self- 
interaction may or may not break parity invariance, and may contribute positively or negatively to the 
energy density, depending on the couplings considered (see [TJ [51 [5J M EH [HI H2] for related literature) . 

In the second part of this paper (Sections [H [5] and [5J we examine cosmological models ruled by a spinor 
field within this theory. We find that the dynamics can be reduced to a closed set of ODEs, representing 
the metric and some of the degrees of freedom contained in the spinor bilinears (Section [4]). We then seek 
solutions to these equations. Solutions exhibiting parity-symmetry (which we label "ambidextrous") are 
particularly simple to integrate. In the presence of torsion, for one sign for a given combination of the 
couplings, we find a bouncing Universe, driven by a torsion-induced phantom phase (Section [5]). In these 
solutions there is a maximal energy density, and the cosmic singularity has been removed. These solutions are 
stable, and range from the very simple to the very complex, as we show for more general, non-ambidextrous 
solutions in Section [51 

In a concluding Section we list open issues to be addressed in future work. We also include two Appendices 
where the notation used in this paper is thoroughly explained. 



2 The theory 

In this paper we will look at the gravitational effect of a specific type of spinorial matter in cosmology. We 
first discuss our choice for the action describing the gravitational field. We shall use a first-order formalism 
for gravity where the gravitational field is described by the spacetime one-forms e 1 = e 1 ^dx^ and oj ij = 
u> IJ fl dx fl . The field e 1 is referred to as the co-tetrad and transforms homogeneously under local SO(l, 3) 
transformations whereas the field uj ij , referred to as the spin-connection, transforms as a gauge field under 
similar transformations. We restrict ourselves to Lagrangians which are generally covariant, locally Lorentz 
invariant, and polynomial in our basic fields and their derivatives. Throughout this paper we shall rely heavily 
on the language of differential forms and we refer the reader to |13| for an introduction to these methods in 
gravitational theory. For compactness of notation we write the wedge product a A b of two differential forms 
simply as ab. The previous requirements on the gravitational Lagrangian restrict the number of possible 
terms considerably. Indeed it may be shown that up to boundary terms, the only 'ingredients' that can 
be used are: co-tetrad e 1 , the curvature of the spin-connection R IJ = duj IJ + uj 1 k uj kj , and the SO(l,3) 
invariant objects €ijkl and r\u = diag(— 1, 1, 1, 1) . Each of these quantities transform homogeneously 
under SO(l,3) transformations and so one can combine these quantities to construct differential forms that 
are Lorentz scalars. We consider the following action: 



S G [e\uj IJ ] = kJ (e IJK L + ^ViKVJL^e I e J R KL (1) 

where k — 1/(32ttG). The first term is the familiar Palatini action, whilst the second term is referred to as 
the Hoist term, with 7 the Immirzi parameter. We do not include a cosmological constant term. The only 
further actions which are functionals only of e 1 and R IJ and polynomial in these fields are boundary terms 
quadratic in R IJ [M] , 

We now consider actions describing spinorial matter. If gravitation is essentially related to local Lorentz 
invariance then this matter will be described by Weyl spinors, i.e. vectors in the left and right handed 
representations of the group £1/(2,(7) [T5]. We may now look to construct the most general action for 
SL(2, C) spinors that produce the familiar spinor equations of motion in flat spacetime. We will use the 
label '(/)' to denote quantities associated with the left handed representation and '(r)' to denote quantities 
associated with the right handed representation. As in the gravitational case, we will look to construct 
actions which are polynomial in fields and locally Lorentz invariant. These actions should contain spacetime 
derivatives of the spinor fields so that they have dynamics. With this in mind, consider the following 
spacetime one-forms constructed from left handed Weyl spinors (f> and right handed Weyl spinors x : 



(2) 
(3) 
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where 



uj ij Cue, 



D^ X = d X ~-uJ 1J KuX 



(4) 
(5) 



The quantities Cu and IZu are the left and right handed generators of SL(2, C) (see Appendix El for 
explicity expressions), whilst the ellipsis in each case denotes terms associated with coupling to different 
gauge fields. For instance, for spinor fields of the standard model of particle physics will contain the 
weak force gauge field whilst will not. Therefore cj) and x mav have additional Yang-Mills indices, 

though we assume that the f operator itself has sufficient additional structure so that the kinetic terms are 
scalars under the relevant Yang-Mills transformations. The one-forms K 1 ^ and K 1 ^' transform as complex 
SO(l, 3) vectors and so in requiring real actions we may consider the following combinations: 



S (0 [^w",e J ] = le IJKL e I e J e K 



S( r )[x,u u ',e Jl 



^UKLe I e J e K 



a (r) {K L{r ^ + K L{r >^j + ib (r) (k l ^ - K L{r >^j 



where the a/, r and 6; >r are real constants. We note by inspection that the following relations hold: 



D{^a I <p) 



Therefore the am and d( r ) terms of the action may be written collectively as follows: 



(G) 
(7) 



(8) 
(9) 



e U KLe I e J e K D{a (l) ^a L (t) + a (r) x^ L x) =3 / (a {l) ^ a L c/> + a {r)X ^ <r L X ) ^ijklT 1 e J e K 

d [eijKL (a (l) tfa L <f> + a {r)X ^ L x) e 1 e J e K ] (10) 



where 



T 1 = De 1 = de 1 + uj 1 je J (11) 

The two-form T 1 is referred to as the spacetime torsion and is assumed to be vanishing in the conventional 
second-order metric formulation of gravity [16]. We see from (jlOl) that the am and at r ) only contribute to 
the equations of motion (i.e. are not described only by a boundary term) when T 1 (x^) is non- vanishing. 
We may additionally consider an action Su [<j>, X> eI ] describing 'potential' terms built from combinations of 
<p and X- 



-~ J U(^x)eijKLe I e J e K e L 



(12) 



Note that the inclusion of factors such as 4! in the above equation and 3! later on are to ensure neatness when 
actions are written in standard form (see Appendix [Bjl . Our total spinor action then is = Sm + + Su, 
i.e. 



£UKLe J e K < ie 1 



1 



+3T 1 (a {l) ^a L ^ + a {r)X ^ L x) - ^Ue I e L 



(13) 
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Additional notational simplification is possible if we describe spinors in terms of Dirac spinors. Following the 
field redefinition X — > y/btn/bt r \X we define the Dirac spinor ^ = (</> Q ,x Q ) an d so the action 5* becomes: 

S 9 [e I ,w IJ ,V] = ^Je I j KL e J e K ^e I (^(^ 1 L D^-D^ 1 L ^+^T I (aV L + (3A L )^ 

~h J U£ i'JKLe I e J e K e L + S int ( X , 4>, ..) (14) 

where "J = ^7° and U is specifically considered a function formed from scalars formed from Dirac spinors 
and the basic objects of spin(l,3) ~ SL(2,C) e.g. *E , \E r or ^7 5 \E r . The covariant derivative is defined 
as D^f = — Ju^ where Ju are the generators of spm(l,3) (see Appendix IA"]) . The action S. m t 
contains interaction terms between (f>, x , and Yang-Mills gauge fields B^. Additionally we have defined the 
following quantities: 

V L = *7 L * = <j)^a L <j) + x^ L X aL = *7 5 7 L * = $ a L <j> - X ^ <j l X 

h - 1 - 1 ( X hl)\ a- 1 ( b (D 



12 12 V w K 'b(r)J 12 v ' ' ' ' o(r) 

The vectors V 1, and A L are referred to, respectively, as the vector and axial current density of the spinor field 
. These quantities will be of particular importance. One may additionally consider terms coupling spinor 
invariants to curvature and torsion, for instance ^hPT^Tr, but we will not consider them in this present 
work. Our combined action will therefore be S = Sc + S*- For ease of comparison with similar actions 
considered in the literature, we present the form of this action in 'standard' tensor notation in Appendix iBl 



3 Equations of motion and first implications 

We now derive the equations of motion. As usual these will be defined by the requirement of stationarity 
of the action S under small variations of the dynamical fields. Varying with respect to the spin-connection 
uj ij we have that: 



2k \e IJKL + -Vi[kVl]j) T'e J + -e IMNP e I e M e N e DP KL A D 
\e [K[MNQ e M e N e [L] (aV Q + pA®) (15) 

Therefore we see that the torsion T 1 is sourced by axial and vector currents A L and V L . For computational 
convenience it will be useful to decompose the spin-connection as follows uj ij — uj ij + C IJ . The quantity 
u> IJ is defined to be the solution to the spin connection when V L — A L — 0, i.e. it is a solution to the 
equation T 1 = 0. Therefore Q IJ depends only upon e 1 and its partial derivatives. By implication the torsion 
may be expressed in terms of the contorsion one-form C IJ = C IJ fi dx fl as follows: T 1 = C 1 ej. Furthermore 
we define the contorsion scalar C IJM eM IJ ,dx ti = C IJ dx^\ after calculation it may be seen that (fT5|) implies 
the following solution for this quantity: 



SS 
S^ 1 



Ctlk 
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in (7 2 + 1 
1 



ytJCg yA D + -(aV D +/3A D ) 
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A[lVt]k + aV[ L T] T ]K + PA [L r] T ] 



K 



(16) 



In summary then we have solved for u IJ in terms of ui IJ and C IJ , which themselves may be expressed 
entirely in terms of ej, V 1 , A 1 and their derivatives. In this sense we have eliminated torsion from the 
theory as the remaining equations of motion may be expressed entirely in terms of variables familiar from 
the second-order formalism. We now find the Einstein equations which follow from varying the action with 
respect to the co-tetrad e 1 : 
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SeJ ~ 



2nej IKL e J R KL - t -e I j KL e J e K (i&j L m - DVj 1 "® 



~e MJ i L e J e M e L + 2k (e JIKL + ^VJKViA C K MP C ML Q e J e p e Q 
+ \eiKNL (~e EL AB AE + \5 L A {aV B + 0A B )j C AB P e N e K e p 

Objects with a'"' above them denote quantities constructed from the 'zero-torsion' spin connection, obtained 
by replacing ui IJ with uj ij wherever they occur within the object. Finally we write down the spinor equation 
of motion which comes from considering variations with respect to 

6 S 

-=-=0 = Aie IJKL e J e K e I 1 L D^ + Qae IJKL e J e K T I 1 L ^ 
ow 

+6 (1 + (3) euBeW^* - %e IJKh e 3 e K e 1 e L (17) 

of 

3.1 The four-fermion interaction 

It is instructive to write the Einstein and Dirac equations in standard tensor notation. By calculation we 
find that: 

-Wg^ (18) 
i 7 MA<* = 6 -Z (19) 
where W is the 'effective potential', incorporating the effects of the non- vanishing contorsion form: 



2 



W = U + 



2(l + 7 2 ) 



1 - (3 2 + ~/3^ A j A 1 - oPViV 1 -2a(p- ~\ AjV 1 



(20) 



(Here we have used the fact that k = l/32irG.) For simplicity we restrict ourselves to the case where 
the potential U depends only on 'F'F. We may use the Dirac equation to simplify the Einstein equations 
somewhat, recasting them in the following form: 

^G, v = -ie i ,(%^-^ 7 ^) + (^ + ^ y ^-2C/)^ (21) 

Equations (|19p and (|2ip are the classical equations of motion for a Dirac spinor field non-minimally coupled 
to gravity. 

Equation (f2"Uf is the central result in the first part of our paper (theoretical set up). It shows that the 
torsion effects predicted by our theory can be recast in the form of 4-fermion self-interactions. This is not 
new (see [6l[TJ[7] for instance), but we have applied this idea to a more general framework. It is useful to 
consider the various limits of the theory. Minimal coupling for the spinor is obtained by removing terms in (J6j) 
and J7} which are pure boundary terms when the torsion vanishes. This amounts to setting am = ai r ) = 0, 
i.e. a = (3 = 0. More generally, when 7 = 0, the self-interaction vanishes, and we recover the torsion-free, 
second order formulation. Indeed setting 7 = in (TT} is equivalent to setting the term multiplying 1/7 to 
zero, i.e. setting the torsion to zero; this is the theory underlying the cosmological models studied in [TTJ. In 
contrast, when 7 —¥ 00 we recover Einstein-Cartan-Sciama-Kibble theory for which the cosmological effect of 
spinorial matter has been considered in some detail jTSJ (HJ . However vector- vector and a parity violating 
axial- vector interactions do appear for more general couplings, as we have shown here. 



() 



In between these extreme cases we find a class of theories parametrized by the Immirzi parameter and 
the non-minimal coupling constants a and (3. The form of W for the case where (3 = has been worked 
out [SI M HD] ! and our general result falls within the results in [TTJ [T2] . The exotic case 7 = ±i sees the 
interaction diverge. This corresponds to setting the (anti)self-dual current to zero in the minimally coupled 
case, and more generally what's inside the bracket in Eq. (|20[). However, our insistance upon real actions 
restricts us to real values for 7. 



3.2 An application: a Classical Spinor Field 

It is expected that the equations of motion (fl9l) and (|21|) should ultimately be regarded as operator equations 
for the spinorial and gravitational quantum fields. We shall simply assume that the classical gravitational 
field that we observe is sourced in (|2"TT) by expectation values of spinor invariants such as ^"3/ and AjA 1 . It was 
noted in [19] that in the first-order formulation of gravity there is an ambiguity in taking expectation values, 
wherein it is arguably more natural to take expectation values in the equation of motion SS/5uj ij = 0. Upon 
solving for the (classical) contorsion, this would yield contributions proportional to (Aj) (A 1 }. If, however, 
one had started from the second-order formalism with four-fermion interaction, it would be expected that 
contributions in the Einstein equations would be of the form (AjA 1 ^. This is an important issue as (Aj) (A 1 ^ 
and (Ai A 1 } will typically not be identical [3T]. We note that our ability to cast our model as a second order 
theory has relied on being able to solve for the contorsion form algebraically; for modest modifications to 
gravity this is no longer possible (see for instance |22j ) and so one may doubt the primacy of the second-order 
formulation of gravity and its implications for the gravitational effect of fermions. 

We circumvent this possible ambiguity by restricting ourselves to spinors called classical spinors. These 
arc defined to be quantum fields described by a spinorial operator ^ = (0, \) & n d assumed to be in a 
state where the expectation value (f{^)) ~ /((4 r )). If this is the case then the above ambiguity in the 
averaging of AjA 1 disappears. Henceforth we will confine ourselves to situations where the approximation 
« is sufficiently good to be regarded as equality, and the field (\I r ) = will be referred to as a classical 
spinor. It should be noted that familiar fields such as those describing quarks and lcptons are not classical 
in this sense on cosmological scales, therefore when we consider a spinor field that is independent from the 
fields of the standard model. If we consider explicit components (in the representation of the Dirac matrices 
given in Appendix A), "I'd = (a, b, c, d), where a, b, c, d are assumed to be complex numbers, we have that: 



(AjA 1 ) = 4(a*c + b*d)(a*c + b*d)* (22) 

{VjV 1 ) = -4(a*c + b*d)(a*c + b*d)* (23) 

(V1A 1 ) = (24) 

= {a*c + b*d) + (a*c + b*d)* (25) 

(vj/7 5 *) = (a*c + b*d) - (a*c + b*d)* (26) 



Therefore, for the classical spinor, a non- vanishing (Ai A 1 } is always spacelike, a non- vanishing (ViV 1 ') is 
always timelike, and they are of equal magnitude. All of the above quantities depend upon a single complex 
number: 



a*c + b*d = i (E + Bi) (27) 



hence we have: 



(A^ 1 ) = - (ViV 1 ) = (E 2 + B 2 ) (28) 
(##) = E (29) 
(* 7 5 *) = iB (30) 

We have explicitly derived these identities in order to motivate the introduction of variables E and B, but 
they could have been obtained more directly with knowledge of the Pauli-Fierz relation: 

(*Q 77 *)(*Q7 / *) = -(*Q*)(*Q*) + (*Q75*)(*Q75*) , (31) 
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where Q G {1,75}. Given the assumption of a classical spinor, the form of the function W simplifies 
considerably: 

W = U{E) +£(E 2 + B 2 ) (32) 

where 

3^G 7 2 (, , 2 , 2 a a2 



2(7 2 + 1) 



1 + a 2 + -13 - (3 Z \ . (33) 



We stress that within the set of couplings considered £ can be positive or negative. Recalling that A 1 must be 
space-like, this means that the contribution to the overall energy due to torsion may be positive or negative 
in our model, a fact that will have far reaching consequences in this paper. 

4 Cosmological equations 

We would like to set up a model based on the FRW metric and a classical spinor field, sourcing torsion. It 
is not immediately obvious that this is possible. Consider the axial current A 1 for a generic classical spinor. 
Since this is space-like, there isn't a frame where its spatial components vanish, and therefore any spinor field 
picks up a preferred direction. However, this does not imply that the metric has to be anisotropic. In fact, 
as Isham and Nelson showed [23], the metric may still be the FRW metric even if the spinor is anisotropic, 
as long as K — 0, i.e. it is the spatially flat FRW metric. Otherwise it is impossible to satisfy the Einstein 
equations, precisely because the A % cannot be made to vanish. We shall therefore assume K = throughout 
this paper. In a sense this is a solution to the flatness problem: the existence of a spinor in a Friedmann 
universe would preclude the existence of spatial curvature. 

The cosmological consequences of a classical spinor have variously been considered in the literature. 
In [T7] the cosmological effect of a classical spinor with potential U(^>^>) and with 7 = (i.e. vanishing 
torsion) up to small perturbations around a spatially flat FRW universe was considered. This analysis has 
subsequently been extended to the Einstein-Cartan minimal coupling case 7 = 00, a — (3 — [2D] , Further 
to this, non-minimal coupling a ^ 0, (3 — has been considered at the level of the cosmological background 
in [23] though we note that our sign for the four-fermion interaction in W does not agree. 

We make the following choice for our co-tetrads: e° = dt, e l = a(t)dx l where indices i, j, k, .. go from 1 to 
3. We first obtain the non- vanishing components of the field ui IJ , the torsion-free spin connection, by solving 
de 1 + Q 1 je 1 — 0. By inspection we have de 1 — HaSfdtdx 1 where H = - and so we have ui' l = He 1 . Given 
our ansatz for the spinor field and geometry, it may be shown that the Dirac equation takes the following 
form: 

and taking the Hermitian conjugate of the above equation we get 



*t + ^*tj 7 o = i^j (35) 

The simple algebraic facts we have just presented are enough to derive a closed set of ordinary differential 
equations for the spinor and gravitational fields. We choose to express the spinor dynamics in terms of its 
quadratic invariants, since these are the observables of the theory (rather than the spinor field itself). Indeed 
it turns out that a complete closed set of equations for our system, assuming the FRW metric with K = 0, 
is formed by: 

E + 3HE = 4£BA° (36) 
B + 3HB = -4£EA° - 2U'(E)A° (37) 
A + 3HA° = 2U'(E)B (38) 

H 2 = ^(U + aE 2 + B 2 ))^p (39) 



5W\ 



t 
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where the prime denotes differentiation with respect to E. This is a minimal set, to which other equations 
could be added. For example: 

V° + 3HV° = (40) 

but this drops out from the gravitational dynamics altogether. Likewise one could write more ODEs ruling 
the dynamics of the remaining variables associated with the bilinears, but they also drop out from the 
cosmological dynamics. For completeness, the second Friedmann equation is given by 

-2--H 2 = 8ttG \U'E - U + £ (E 2 + B 2 )] = SitGp (41) 
a L 

This equation can be derived from (|39[) and conservation equation 

p + 3H{p + p)=0 (42) 

which in turn follows from the (|36l) . (|3"T| and (|3"5|) . This is, however, not true at turn-around points (where 
H = 0), because the conservation equation becomes degenerate (0 = 0). The second Friedmann equation is 
then needed, and failure to take this fact into account might lead one to mistake a bounce for a static or 
loitering Universe. Finally we note that equations ((36)) . (|37|) and (138)) provide a first integral: 

a 6 [E 2 + B 2 + (A ) 2 )] =M 2 (43) 

This equation was first noted in [2~i] . 



5 Ambidextrous solutions 

We can immediately identify a number of possible effects of the spinor that persist even if we add an extra 
fluid to the system. Single-chirality Weyl spinors, <fr or x, for example, are not very interesting in cosmology. 
They produce vanishing currents A 1 and V 1 , with E = B = 0. Such spinors therefore have no effect on 
the Friedmann equations, and if there are no other matter components in the Universe, they simply lead 
to Minkowski space-time. Single-entry spinors in the Dirac representation, on the other hand, map into 
solutions with \(f)\ 2 = \x\ 2 , i.e. solutions without parity violation, where left and right spinors have the 
same probability. These "ambidextrous" solutions do have an effect on the Friedmann equations. They are 
particularly simple to integrate because they display 

A = (\<f>\ 2 -\x\ 2 )=0 (44) 

which at least when U' ^ implies B = (see Eqs. (l37l) and (|38|l). Then Eqn. (|36f provides the first 
integral: 

or 

which is nothing but (l4"3l when we set B = A° =0. 

As we see, we must distinguish between parity violation in the solutions to the theory, and in the theory 
itself (i.e. in its action), here only present if (3 / I/7 (cf. the last term in Eqn. (|20l0 . Regardless of the 
parameters of the theory we see that solutions with maximal parity violation (Weyl spinors) have no effect 
in cosmology, with or without torsion, and with or without parity violation in the actual theory. Parity 
invariant or ambidextrous solutions, in contrast, have an effect particularly simple to analyze, which we shall 
now do. In between the two extremes we find intermediate solutions, harder to work out, which we will do 
numerically in Section [5] 

5.1 Solutions without torsion 

If 7 = 0, it follows that £ = 0, and we obtain the well-known case of the torsion- free theory. As pointed out 
in |17) , classical spinors are remarkable (and bypass a number of theorems valid for scalar fields) in that any 
equation of state can be produced by appropriately designing the potential. For example, if: 

U = ((m) n = (E n (46) 
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then, since p = U and E oc 1/a 3 , we have p oc l/a 3 ™. 
state w — p/p) p oc l/a 3 ^ 1+w ^ we can read off 



Since (|42|) implies that (for a constant equation of 



Wq = n — 1 



(47) 



without any further calculation. The standard results for a(t) follow (e.g. a oc 

£2/(31+^ Jf y, > 

The n = 1 case corresponds to the massive Dirac field, whereas n = 4/3 corresponds to the Gursey model 
[25] , In the absence of self-interactions, the former leads to a dust model, whereas the latter behaves like a 
radiation dominated Universe. Inflation can only be precisely obtained with a flat potential (n — 0), with 
the problems discussed in [I?] . 

5.2 Torsion driven bouncing solution 




Figure 1: The bounce for a typical ambidextrous solution with £ < 0. In this case U = mE, so that the 
contracting phase exhibits a oc (— t) 2 / 3, and the expanding phase a oc < 2 / 3 . 



If 7 7^ in general we have a bouncing solution. For simplicity let us consider the n = 1 mass potential, 
U = mE, but the results in this and the next sections generalize to more complex potentials. The case £ > 
will lead to less interesting solutions, to be reviewed in Section IfTSl If £ < 0, the self-interactions resulting 
from torsion have the effect of dramatically reshaping the mass potential, converting the typical mass bowl 
into an upside-down Mexican hat potential: 



W = mE- \£,\E 2 . 



(48) 



Since p = W > (as implied by the first Friedman equation, Eq. ([31?))') it is then not difficult to predict a 
bounce when p — 0. The only alternative would be a static universe, but since p^O when p = this is not 
realized. 

We have plotted a(t) in Fig. [TJ Asymptotically (large |i|), the Universe contracts like a oc (— i) 2 / 3 and 
expands like a oc i 2 / 3 , typical of dust. In between there is a phase where torsion dominates causing a bounce. 
This occurs when: 

(49) 
(50) 



E — E — — — . 



We may write 



P = P 



Po 



where p = mE represents the energy density before torsion effects are added, and 



Po = 



(51) 
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Then we see that the bounce occurs when p = and p = pq. By studying further the function p{po), we see 
that it has a maximum at 

Pmax = ■ (52) 

As the universe contracts the density increases like 1/i 2 , as usual, but then torsion kicks in. This maximal 
density is reached, and then, as the universe compresses further, the density decreases until it reaches zero 
and a bounce occurs at a finite a. After the bounce the density at first increases with expansion, until the 
same maximum p max is reached again. After that it starts to decrease with expansion, eventually according 
to the usual p oc 1/i 2 . The density never diverges and a Big Bang singularity is avoided. We have plotted 
this behaviour in Fig. [2] 




Figure 2: The values of p, displayed in units of po, and w during an ambidextrous bounce. As we can see 
there is a period of phantom behaviour, where the density decreases as the universe contracts, then increases 
as it expands. 

Increasing density with expansion (or decreasing density with contraction) is a hallmark of phantom 
matter [55]. And indeed by plotting the equation of state w = p/p (see Fig. [5]) we see that this does cross 
w = — 1 as the maximum density is reached. Torsion induces a phantom period around the bounce and in 
fact w becomes infinitely negative at the bounce, since p < is finite while p = 0. When torsion becomes 
sub-dominant w goes to zero, as predicted in the previous sub-section. 

6 Parity violating perturbations 

The issue arises as to whether the bouncing solution presented in the previous Section is stable, when strict 
parity invariance is broken, with B and ^4° turned on. We will study the matter in this Section, first 
turning on £?-type perturbations at the bounce, then A perturbations, and then both. In all cases the 
solution is stable, in the sense that we find small variations in the bounce commensurate with the size of the 
perturbations induced. What is more important, we find that even when the perturbations are very large 
the overall picture does not change much in the first two cases (pure B and ^4° type perturbations), the 
solutions simply displaying (large) oscillations around the basic bouncing solution. However, when B and A 
perturbations are allowed free rein at the bounce they introduce a very interesting qualitative novelty should 
they be large enough: whilst the bounce is still present there is an asymmetry between the contracting and 
expanding phase. 

6.1 B-type and y4°-type perturbations 

As just stated, we do not see any instability in the bouncing solution, and small additions of B and A 
lead to perturbations of the same order. What is interesting is that even if the perturbations are very large 
(order 1 and higher) the overall picture does not qualitatively change, as long as one of B and ^4° are zero 
at the bounce. If B = at the bounce we call this an A°-type perturbation and vice- versa. In Fig. [3] we 
plot the effect of a large -B-type perturbation on p and w. The overall picture is essentially the same, with 
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Figure 3: The values of p and w with a large 5-type perturbation at the bounce. Even though the a(t) 
profile doesn't change much we find interesting oscillatory behaviour superposed on the basic, ambidextrous 
picture. 



oscillatory behaviour superposed on the ambidextrous picture. This is due to the oscillatory nature of the 
variables B and ^4° (see Fig. |H top). Obviously the value for the maximal density is now different, and most 
easily determined numerically, but the basic picture remains. The bounce occurs when 



E = 



—m ± y/m 2 - 4£ 2 5 2 
2£ 



(53) 



The picture is similar for an A°-type perturbation. In Fig. 0] we plot the behaviour of B and A for .B-type 
and A-type large perturbations, for comparison. The two variables are generally out of phase, and the 
modes considered here correspond to one of them having a node at the bounce. We don't plot the a(t) 
profile because this is basically indistinguishable form the unperturbed case described in Fig. [TJ 



B AO 

1.0 r 



B AO 

1.0 r 




Figure 4: A plot of B (blue) and A (red) for B-type and A-type large perturbations. 



Another way to represent these results is to map the dynamics onto the plane spanned by 

a 3 B 



(54) 



M 
a3A ° 

V = ~W (55) 

Given Eq. (|43[) . the system is constrained to remain inside the unit circle, with the distance to the boundary 
providing a measure of a 3 E. In this plot, the origin corresponds to the exactly ambidextrous solution 
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presented presented in the last Section. The {x, y} trajectories away from the origin represent the parity 
violating perturbations. In Fig. [5] and |5] we show these trajectories for i?-type and A-type perturbations, 
respectively. The outer trajectories correspond to the rather large perturbations used in the previous plots. 
The inner trajectories correspond to smaller and smaller perturbations in the initial conditions. As we see 
the ambidextrous solution is stable; but more interestingly even very large perturbations arrange themselves 
as oscillations around this solution, which therefore seems a generic feature. 



Figure 5: The {x, y} trajectories of B-type perturbations (trajectories constrained to remain inside the unit 
circle, depicted). The outer trajectory corresponds to the large perturbations used in the previous plots. The 
center represents the exactly ambidextrous solution presented in the last Section. The intermediate trajec- 
tories correspond to smaller and smaller perturbations in the initial conditions. As we see the ambidextrous 
solution is stable, but furthermore, even very large perturbations arrange themselves as oscillations around 
this solution. 



6.2 Generic perturbations 

The only qualitative novelty appears if we mix the two types of perturbations, and let their amplitude be 
very large. If we set up the perturbations away from the bounce in general it does not happen that one of 
B or A vanish at the bounce. Then, it still happens that the ambidextrous solution is stable against small 
perturbations, and that even very large perturbations consist of oscillations around this solution. However a 
novelty appears: there appears an asymmetry between the contracting and expanding phase, as we illustrate 
in Fig. [71 For example, the maximal energy density reached is different in the two phases, as is the defining 
"constant" a 3 p associated with the two dust Universes. Depending on the initial conditions chosen, this may 
be smaller or larger in the expanding phase. We also depict the evolution in an {x,y} plot in Fig. [5] The 
evolution is still an oscillation around the ambidextrous solution (origin) but now the amplitude on cither 
side of the bounce is different. 

6.3 Other solutions 

There are a number of other solutions which are not particularly interesting, but which we list here for 
completeness. 



a 



3 A 



M 
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Should £ > 0, then there isn't a bounce and the singularity is not avoided. Indeed as we go back in time, 
the singularity is reached faster because the torsion adds to the pressure. The dust phase described above 
is then preceded by a period of kination (to borrow terminology from scalar field cosmologies), i.e. a period 
with p oc 1/a 6 and a oc i 1 / 3 . This is true for ambidextrous solutions, but does not qualitatively change for 
more general solutions. If the theory is massless there is only a period of kination, and this happens for the 
more general case with non- vanishing B. In that case B oc 1/a 3 and it does not affect this conclusion. It 
could also be that U = Uq is a constant, in which case the period of kination with E and B decaying like 
1/a 3 is followed by a de Sitter phase. None of this is very surprising or interesting. The strength of this 
paper is in the £ < solutions. 



7 Conclusions 

We conclude with an appraisal of what we have discovered and what remains to be done. We found that the 
overall picture of cosmologies driven by spin hinges on the sign of one combination of parameters, £ as defined 
in Eq. (1331) . If £ > 0, as is the case of Einstein-Cartan theory with minimal coupling, nothing dramatically 
new happens. For example, for a mass potential (or any potential leading to w < —1/3) the only novelty 
is that an early period of kination precipitates the onset of a singularity as we go back in time. But if 
£ < the singularity is avoided and even for a simple mass potential a bounce is generic, even against very 
large perturbations of the basic solution. This double picture is closely related to the attractive/repulsive 
nature of the 4-fermion interaction. As is well known, in the Einstein-Cartan theory with minimal coupling 
this interaction is attractive (at least for classical spinors), with the result that singularities are easier to 
form than in the torsion-free theory [27j . By identifying a non-minimal coupling leading to £ < we have 
reversed the situation, rendering the spinor self-interaction repulsive, and opening up the doors to singularity 
avoidance. It is curious to note that such scenarios are closely related to the presence of parity-violating 
terms in the action. 

Several questions can be raised. We presented an extensive set of solutions, displaying simple and 
complicated bounces which may be symmetric or not. But have we found all possible solutions? It is 
tempting to speculate on the existence of a de Sitter/inflationary fixed point. Indeed the effects of torsion 
when £ < and fi / are qualitatively similar to those of the C-field in steady-state cosmology: a sea 
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Figure 7: For general large perturbations (where neither B nor ^4° have a node at the bounce) we observe 
an asymmetry between the contracting and expanding phase, with the maximal density different in the two 
phases. 



a 3 AO 
M 




Figure 8: The {x,y} diagram for the evolution of the system under a generic large parity violating pertur- 
bation. We see that the two circles, corresponding to the period when torsion has died down, differ for the 
contracting and expanding phase. The evolution during the bounce is also more complex and asymmetric. 



of negative energy (here represented by torsion) which interacts with a positive energy component (when 
B^O and A ^ 0). Such situations often lead to inflation, i.e. a sustained period of accelerated expansion. 
We have been unable to find such a solution. The reason is probably that in our case the same field has both 
a positive and a negative energy component, so the situation is not quite the same. Nonetheless we defer to 
a future paper a more complete analysis, based on phase space portraits of autonomous dynamical systems. 
We stress that there is accelerated expansion at the bounce (a > 0) but this is not inflation, which is a 
sustained period of accelerated expansion (this seems to have been missed in [IE])- The possibility cannot 
be dismissed, however, that a long inflationary transient is present somewhere in the phase space. 
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Is the bounce we have found stable against anisotropy domination? This is a valid concern because even 
for small perturbations the shear tensor cry contributes to the Friedmann equation as a 2 — <7y(T* J /2, and 
a oc 1/a 3 . Therefore, unless the equation of state is super-stiff (w > 1) the anisotropy tends to dominate 
during the contraction, leading to a mixmaster phase or even a singularity, rather than a smooth bounce. 
In our case we are in a borderline situation, since if torsion and shear do not interact then they both 
scale as 1/a 6 . Therefore the solution we have presented is certainly stable against small, but not large 
shear perturbations. However the situation may be more subtle, and we defer further analysis to a future 
publication. As we noted before, a spinor field already is anisotropic; however this need not be reflected in 
the metric (at least if K = 0). 

Even ignoring the issue of anisotropy, the obvious next step is to work out the fluctuations in this type of 
model, finding the amplitude, spectral index and tensor / scalar ratio as a function of the free parameters of the 
model. This has been examined in the past for spinor-driven cosmologies in the context of inflation [171 120] 
(which may be achieved by making U very flat). As the work of [17] shows, spinors and scalar fields are 
very different in this respect: for the same background kinematics one gets a spectral index n$ = 4 for the 
former where for the latter n$ — 1 is found. For this reason the status of our model regarding fluctuations 
is far from obvious. It is interesting to point out that for scalar fields a w = (dust-like) bounce does have 
a scale-invariant mode [29l [30] . However these models suffer from fine-tuning problems related to the fact 
that the spectrum of the curvature and potential fluctuations is not the same. In future work we hope to 
examine how this might change if the bounce were to be driven by a spinor field. 

The major concern remains as to what might be the physical basis for spinor models (but this criticism 
could be levelled at most early Universe models, including those based on scalar fields). It obviously would 
be more conservative to take the "spin-fluid" approach, such as that pioneered by Weyssenhoff (see for exam- 
ple [311 132 ESI US] for some very interesting cosmological work based on this approach^) , or by considering 
the cosmological effect of thcrmalizcd fcrmionic matter jTHl [XS] ■ However, it may be argued that it is also 
legitimate to consider a spinor field as envisaged here (see Appendix of [17j . for example). In addition we 
note that the theory can be phrased wholly in terms of the bilinear invariants (such as the ODEs presented 
in Section [?]). With this remark in mind in future work we hope to refine the argument in [17j . Notice that 
the sign issues presented in Section 13.21 change dramatically depending on whether the underlying field can 
be considered to be classical (for example ty^iff is positive definite classically but of course need not 

be so in the quantum theory). We close by pointing out that the model presented here has equations very 
similar to those found in loop quantum cosmology and the brane-world scenario, but only when B = 0. If 
B =/= we seem to be generalizing the dynamics in those models. 
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A Conventions 

In this section we describe the conventions used in this paper for various quantities. Often this information 
is omitted in papers about the role of torsion in cosmology, making comparison of results difficult. 

We use distinct symbols hijkl and £^ aS . The object cjjkl is a spacetime scalar antisymmetric in all 
indices and with €0123 = 1; the object is invariant under local 50(1,3) transformations. The object e^ vaS 
is a spacetime density antisymmetric in all indices and with e 0123 = 1; the object is numerically invariant 
under local coordinate transformations. The determinant e of the co-tetrad e 1 = e^dx^ is subsequently 
defined as follows: 

e = ^e IJKL e^elele?et (56) 

1 Note though that there appears to be an essential friction between the Weyssenhoff spin fluid and any underlying theory 
based on minimally coupled fcrmions. The reason for this is that the object Cijk does not possess the same symmetries in 
the two cases 1181 
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and hence 



e^e^ef e L a = ee lJKL = -ee 1JKL (57) 

Furthermore we use the convention that rju = diag(— 1, 1, 1, 1) which implies that the spacetime metric 
Qtiv — Vu e ji e t nas mostly positive signature. In using spinors in this paper we will use the Weyl/chiral 
representation i.e. 



* - I ? ) < 58 > 

where a and a' are indices of the left and right handed representation of SL(2,C) respectively [T5] - We 
choose the following convention for the gamma matrices r ) 1 : 

7 fcV° ' 7 I 1 



a 1 = (l,-a% 
a 1 = (1,0 

where the a 1 are the Pauli sigma matrices. The spin(l,3) generators J IJ take the following form: 

J U = -{[l 1 ,^] (59) 
We now detail our curvature conventions. Our basic object representing curvature is the two-form R IJ : 

R IJ = duo 1 ' 1 +u i k lo kj (60) 
We can define the 'orthonormal components' R IJ ' KL via the following relation: 

R IJ = ^R IJ '^dx^ A dx u = -R IJ KL e^e^dx" A dx u (61) 

= \R IJ KL e K e L (62) 

We define the Ricci tensor IZ^ as follows: 

Tl^ = e Ilx e J v R IK JK (63) 
Consequently the Ricci scalar TZ = R£ is defined as: 

K = R JK JK (64) 

and we define the Einstein tensor as: 

G,a, = H^u - \llg liv (65) 
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B Actions In Standard Notation 

We now use the results of the previous section to write our actions in conventional notation. This will make 
it easier to make contact with antecedent results in the literature. We first begin with the gravitational 
action: 

S G = kJ (e IJKL + ^tiikVJl) e I e J R KL 

= I J (euKL+ 2 -ViKVJi)R KL M N e I e J e^e N 



[eiJKL + ^VlKVJL^j R KL MN^iefe^ dx» A dx v A dx 5 A dx° 
(eiJKL + ^ViKVJL^ R KL MN ey y efe^e^d A x 
cijkl + -VikVjl) R KL MN e IJMN ed 4 x 



K 

2 

u{ eIJKL+2 J 



MN pKL ,4 
KL I n MN ea x 



'■' / I 2R MN MN ~t~ ~^K[j MN R^^MN ^ ' 
1 



/ j 2R + ^e KL MN R KL MN )< l r r 



Clearly then we have k = l/327rG. Next we write the spinor action in a more familiar form. We first 
consider the the action in the limit of zero torsion T 1 = 0: 

Sa, iT i= = y ^e IJK Le I e J e K ^ L D^ - —e IJKL e I e J e K e L - h.c 

= h j [^JKLeieief^ L D^ - |e /JA r L e^e^f e^j e^d 4 x - h.c 



3! 

^ / (y.iKLe IJKM el I ^ L D a ^-^ I . !K L: ,IKL )-/ , --/,, 



2< Jf e^^Da* - U ) ed^x - h.c 



e1^j L D^ - U ) ed 4 x - h.c 



Finally we turn to the non-minimal coupling terms of the left hand side of equation (|10[) . In terms of our 
Dirac spinor <J/ (and recalling the redefinition of \ following (|10l0 this action, S nm can be written: 

S nm = ^ J e IJKL c I e J e K (a^ L DV + 0j 5 ^ L D^>) + h.c 
= ~ J el (a^ L D^ + /3*7 5 7 L L> M *) ed 4 x + h.c 
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